Human gene regulatory evolution is driven by the divergence of regulatory element function in both cis and trans.

Gene regulatory divergence between species can result from cis-acting local changes to regulatory element DNA sequences or global trans-acting changes to the regulatory environment. Understanding how these mechanisms drive regulatory evolution has been limited by challenges in identifying trans-acting changes. We present a comprehensive approach to directly identify cis- and trans-divergent regulatory elements between human and rhesus macaque lymphoblastoid cells using assay for transposase-accessible chromatin coupled to self-transcribing active regulatory region (ATAC-STARR) sequencing. In addition to thousands of cis changes, we discover an unexpected number (∼10,000) of trans changes and show that cis and trans elements exhibit distinct patterns of sequence divergence and function. We further identify differentially expressed transcription factors that underlie ∼37% of trans differences and trace how cis changes can produce cascades of trans changes. Overall, we find that most divergent elements (67%) experienced changes in both cis and trans, revealing a substantial role for trans divergence-alone and together with cis changes-in regulatory differences between species.

Cross-tissue patterns of DNA hypomethylation reveal genetically distinct histories of cell development.

BACKGROUND: Establishment of DNA methylation (DNAme) patterns is essential for balanced multi-lineage cellular differentiation, but exactly how these patterns drive cellular phenotypes is unclear. While > 80% of CpG sites are stably methylated, tens of thousands of discrete CpG loci form hypomethylated regions (HMRs). Because they lack DNAme, HMRs are considered transcriptionally permissive, but not all HMRs actively regulate genes. Unlike promoter HMRs, a subset of non-coding HMRs is cell type-specific and enriched for tissue-specific gene regulatory functions. Our data further argues not only that HMR establishment is an important step in enforcing cell identity, but also that cross-cell type and spatial HMR patterns are functionally informative of gene regulation. RESULTS: To understand the significance of non-coding HMRs, we systematically dissected HMR patterns across diverse human cell types and developmental timepoints, including embryonic, fetal, and adult tissues. Unsupervised clustering of 126,104 distinct HMRs revealed that levels of HMR specificity reflects a developmental hierarchy supported by enrichment of stage-specific transcription factors and gene ontologies. Using a pseudo-time course of development from embryonic stem cells to adult stem and mature hematopoietic cells, we find that most HMRs observed in differentiated cells (~ 60%) are established at early developmental stages and accumulate as development progresses. HMRs that arise during differentiation frequently (~ 35%) establish near existing HMRs (≤ 6 kb away), leading to the formation of HMR clusters associated with stronger enhancer activity. Using SNP-based partitioned heritability from GWAS summary statistics across diverse traits and clinical lab values, we discovered that genetic contribution to trait heritability is enriched within HMRs. Moreover, the contribution of heritability to cell-relevant traits increases with both increasing HMR specificity and HMR clustering, supporting the role of distinct HMR subsets in regulating normal cell function. CONCLUSIONS: Our results demonstrate that the entire HMR repertoire within a cell-type, rather than just the cell type-specific HMRs, stores information that is key to understanding and predicting cellular phenotypes. Ultimately, these data provide novel insights into how DNA hypo-methylation provides genetically distinct historical records of a cell's journey through development, highlighting HMRs as functionally distinct from other epigenomic annotations.

A role for N6-methyldeoxyadenosine in C. elegans mitochondrial genome regulation.

Epigenetic modifications provide powerful means for transmitting information from parent to progeny. As a maternally inherited genome that encodes essential components of the electron transport chain, the mitochondrial genome (mtDNA) is ideally positioned to serve as a conduit for the transgenerational transmission of metabolic information. Here, we provide evidence that mtDNA of C. elegans contains the epigenetic mark N6-methyldeoxyadenosine (6mA). Bioinformatic analysis of SMRT sequencing data and methylated DNA IP sequencing data reveal that C. elegans mtDNA is methylated at high levels in a site-specific manner. We further confirmed that mtDNA contains 6mA by leveraging highly specific anti-6mA antibodies. Additionally, we find that mtDNA methylation is dynamically regulated in response to antimycin, a mitochondrial stressor. Further, 6mA is increased in nmad-1 mutants and is accompanied by a significant decrease in mtDNA copy number. Our discovery paves the way for future studies to investigate the regulation and inheritance of mitochondrial epigenetics.

ATAC-STARR-seq reveals transcription factor-bound activators and silencers across the chromatin accessible human genome.

Massively parallel reporter assays (MPRAs) test the capacity of putative gene regulatory elements to drive transcription on a genome-wide scale. Most gene regulatory activity occurs within accessible chromatin, and recently described methods have combined assays that capture these regions, such as assay for transposase-accessible chromatin using sequencing (ATAC-seq), with self-transcribing active regulatory region sequencing (STARR-seq) to selectively assay the regulatory potential of accessible DNA (ATAC-STARR-seq). Here, we report an integrated approach that quantifies activating and silencing regulatory activity, chromatin accessibility, and transcription factor (TF) occupancy with one assay using ATAC-STARR-seq. Our strategy, including important updates to the ATAC-STARR-seq assay and workflow, enabled high-resolution testing of ~50 million unique DNA fragments tiling ~101,000 accessible chromatin regions in human lymphoblastoid cells. We discovered that 30% of all accessible regions contain an activator, a silencer or both. Although few MPRA studies have explored silencing activity, we demonstrate silencers occur at similar frequencies to activators, and they represent a distinct functional group enriched for unique TF motifs and repressive histone modifications. We further show that Tn5 cut-site frequencies are retained in the ATAC-STARR plasmid library compared to standard ATAC-seq, enabling TF occupancy to be ascertained from ATAC-STARR data. With this approach, we found that activators and silencers cluster by distinct TF footprint combinations and these groups of activity represent different gene regulatory networks of immune cell function. Altogether, these data highlight the multi-layered capabilities of ATAC-STARR-seq to comprehensively investigate the regulatory landscape of the human genome all from a single DNA fragment source.

Endogenous localization of TOP-2 in C. elegans using a C-terminal GFP-tag.

To investigate the dynamic localization of Topoisomerase II in live C. elegans we have generated a C-terminally GFP-tagged version of TOP-2 at the endogenous locus. We found that TOP-2::GFP localizes in a similar pattern to the previously published TOP-2::3XFLAG strain and does not disrupt the meiotic chromosome segregation functions of this enzyme.

Impact of WIN site inhibitor on the WDR5 interactome.

The chromatin-associated protein WDR5 is a promising pharmacological target in cancer, with most drug discovery efforts directed against an arginine-binding cavity in WDR5 called the WIN site. Despite a clear expectation that WIN site inhibitors will alter the repertoire of WDR5 interaction partners, their impact on the WDR5 interactome remains unknown. Here, we use quantitative proteomics to delineate how the WDR5 interactome is changed by WIN site inhibition. We show that the WIN site inhibitor alters the interaction of WDR5 with dozens of proteins, including those linked to phosphatidylinositol 3-kinase (PI3K) signaling. As proof of concept, we demonstrate that the master kinase PDPK1 is a bona fide high-affinity WIN site binding protein that engages WDR5 to modulate transcription of genes expressed in the G2 phase of the cell cycle. This dataset expands our understanding of WDR5 and serves as a resource for deciphering the action of WIN site inhibitors.

ATAC-Me Captures Prolonged DNA Methylation of Dynamic Chromatin Accessibility Loci during Cell Fate Transitions.

DNA methylation of enhancers is dynamic, cell-type specific, and vital for cell fate progression. However, current models inadequately define its role within the hierarchy of gene regulation. Analysis of independent datasets shows an unanticipated overlap between DNA methylation and chromatin accessibility at enhancers of steady-state stem cells, suggesting that these two opposing features might exist concurrently. To define their temporal relationship, we developed ATAC-Me, which probes accessibility and methylation from single DNA library preparations. We identified waves of accessibility occurring rapidly across thousands of myeloid enhancers in a monocyte-to-macrophage cell fate model. Prolonged methylation states were observed at a majority of these sites, while transcription of nearby genes tracked closely with accessibility. ATAC-Me uncovers a significant disconnect between chromatin accessibility, DNA methylation status, and gene activity. This unexpected observation highlights the value of ATAC-Me in constructing precise molecular timelines for understanding the role of DNA methylation in gene regulation.

Localized TWIST1 and TWIST2 basic domain substitutions cause four distinct human diseases that can be modeled in Caenorhabditis elegans.

Twist transcription factors, members of the basic helix-loop-helix family, play crucial roles in mesoderm development in all animals. Humans have two paralogous genes, TWIST1 and TWIST2, and mutations in each gene have been identified in specific craniofacial disorders. Here, we describe a new clinical entity, Sweeney-Cox syndrome, associated with distinct de novo amino acid substitutions (p.Glu117Val and p.Glu117Gly) at a highly conserved glutamic acid residue located in the basic DNA binding domain of TWIST1, in two subjects with frontonasal dysplasia and additional malformations. Although about one hundred different TWIST1 mutations have been reported in patients with the dominant haploinsufficiency Saethre-Chotzen syndrome (typically associated with craniosynostosis), substitutions uniquely affecting the Glu117 codon were not observed previously. Recently, subjects with Barber-Say and Ablepharon-Macrostomia syndromes were found to harbor heterozygous missense substitutions in the paralogous glutamic acid residue in TWIST2 (p.Glu75Ala, p.Glu75Gln and p.Glu75Lys). To study systematically the effects of these substitutions in individual cells of the developing mesoderm, we engineered all five disease-associated alleles into the equivalent Glu29 residue encoded by hlh-8, the single Twist homolog present in Caenorhabditis elegans. This allelic series revealed that different substitutions exhibit graded severity, in terms of both gene expression and cellular phenotype, which we incorporate into a model explaining the various human disease phenotypes. The genetic analysis favors a predominantly dominant-negative mechanism for the action of amino acid substitutions at this highly conserved glutamic acid residue and illustrates the value of systematic mutagenesis of C. elegans for focused investigation of human disease processes.

The Identification of a Novel Mutant Allele of topoisomerase II in Caenorhabditis elegans Reveals a Unique Role in Chromosome Segregation During Spermatogenesis.

Topoisomerase II alleviates DNA entanglements that are generated during mitotic DNA replication, transcription, and sister chromatid separation. In contrast to mitosis, meiosis has two rounds of chromosome segregation following one round of DNA replication. In meiosis II, sister chromatids segregate from each other, similar to mitosis. Meiosis I, on the other hand, segregates homologs, which requires pairing, synapsis, and recombination. The exact role that topoisomerase II plays during meiosis is unknown. In a screen reexamining Caenorhabditis elegans legacy mutants isolated 30 years ago, we identified a novel allele of the gene encoding topoisomerase II, top-2(it7). In this study, we demonstrate that top-2(it7) males produce dead embryos, even when fertilizing wild-type oocytes. Characterization of early embryonic events indicates that fertilization is successful and sperm components are transmitted to the embryo. However, sperm chromatin is not detected in these fertilized embryos. Examination of top-2(it7) spermatogenic germ lines reveals that the sperm DNA fails to segregate properly during anaphase I of meiosis, resulting in anucleate sperm. top-2(it7) chromosome-segregation defects observed during anaphase I are not due to residual entanglements incurred during meiotic DNA replication and are not dependent on SPO-11-induced double-strand DNA breaks. Finally, we show that TOP-2 associates with chromosomes in meiotic prophase and that chromosome association is disrupted in the germ lines of top-2(it7) mutants.

Discovery of two GLP-1/Notch target genes that account for the role of GLP-1/Notch signaling in stem cell maintenance.

A stem cell's immediate microenvironment creates an essential "niche" to maintain stem cell self-renewal. Many niches and their intercellular signaling pathways are known, but for the most part, the key downstream targets of niche signaling remain elusive. Here, we report the discovery of two GLP-1/Notch target genes, lst-1 (lateral signaling target) and sygl-1 (synthetic Glp), that function redundantly to maintain germ-line stem cells (GSCs) in the nematode Caenorhabditis elegans. Whereas lst-1 and sygl-1 single mutants appear normal, lst-1 sygl-1 double mutants are phenotypically indistinguishable from glp-1/Notch mutants. Multiple lines of evidence demonstrate that GLP-1/Notch signaling activates lst-1 and sygl-1 expression in GSCs within the niche. Therefore, these two genes fully account for the role of GLP-1/Notch signaling in GSC maintenance. Importantly, lst-1 and sygl-1 are not required for GLP-1/Notch signaling per se. We conclude that lst-1 and sygl-1 forge a critical link between Notch signaling and GSC maintenance.